One-neutron removal reactions on light neutron-rich nuclei 
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A study of high energy (43-68 MeV/nucleon) one- neutron removal reactions on a range of neutron- 
rich psd-shell nuclei (Z = 5-9, A = 12-25) has been undertaken. The inclusive longitudinal and 
transverse momentum distributions for the core fragments, together with the cross sections have 
been measured for breakup on a carbon target. Momentum distributions for reactions on tantalum 
were also measured for a subset of nuclei. An extended version of the Glauber model incorporating 
second order noneikonal corrections to the JLM parametrisation of the optical potential has been 
used to describe the nuclear breakup, whilst the Coulomb dissociation is treated within first order 
perturbation theory. The projectile structure has been taken into account via shell model calcula- 
tions employing the psd-interaction of Warburton and Brown. Both the longitudinal and transverse 
momentum distributions, together with the integrated cross sections were well reproduced by these 
calculations and spin-parity assignments are thus proposed for 15 B, 17 C, 19 ~ 21 N, 21 ' 23 , 23 ~ 25 F. In 
addition to the large spectroscopic amplitudes for the ^2sx/ 2 intruder configuration in the N = 9 
isotones, 14 B and 15 C, significant u2s^ 2 admixtures appear to occur in the ground state of the neigh- 
bouring N = 10 nuclei 15 B and 16 C. Similarly, crossing the N=14 subshell, the occupation of the 
z/2si/2 orbital is observed for 23 0, 24 > 25 F. Analysis of the longitudinal and transverse momentum 
distributions reveals that both carry spectroscopic information, often of a complementary nature. 
The general utility of high energy nucleon removal reactions as a spectroscopic tool is also examined. 

PACS number(s): 24.10.-i, 25.60.-t, 25.60.Dz, 25.60.Gc 
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I. INTRODUCTION 



High energy heavy- ion projectile fragmentation has been investigated now for some 25 years [1,2]. The initial 
emphasis centred on reactions employing stable beams at relativistic energies [3,4] and the observed momentum 
distributions (typically Gaussian in form with FWHM <~ 200 MeV/c) were interpreted within a statistical description 
as reflecting the Fermi momentum of the removed nucleons [5] . Further refinements lead to more sophisticated models 
incorporating the peripheral nature of the reaction process [6,7]. Beyond relatively simple considerations, such as the 
surface cluster model [7], the projectile structure played no role in the fragment distributions. 

More recently, the investigation of fragmentation reactions using radioactive beams has lead to the recognition 
of narrow fragment momentum distributions (FWHM <~ 50 MeV/c) as a signature of the extended valence nucleon 
density distributions in halo nuclei [8]. Originally these momentum distributions were assumed, in terms of the 
transparent limit of the Serber model [9], to be a direct mapping via the Fourier transform, of the valence nuclcon(s) 
wave function. Detailed investigation of the role played by the reaction mechanisms, particularly in the case of the 
single-neutron halo nucleus 11 Be [10,11] demonstrated that such an interpretation was not generally applicable. The 
longitudinal core fragment momentum distributions have been suggested [12-14] to be less influenced by such effects 
and as such to represent a cleaner probe. As foreseen in the original work of Serber [9] and more explicitly in the case 
of heavy ions by Hufner and Nemes [6] , the requirement of core survival drives few nucleon removal to probe essentially 
only that part of the valence nucleon(s) wave function residing outside the core [17,18]. Consequently, in the spirit of 
the spectator-core model of Hussein and McVoy [19], various Glauber- type approaches to modelling the dissociation 
of energetic beams of nuclei far from stability have been developed [15-18,20-25]. The essential results are that the 
momentum distributions, as first recognized by Bonaccorso and Brink [26] and Sagawa and Yazaki [15], reflect the 
angular momentum (I) of the removed nucleon, whilst the corresponding cross section may provide a measure of the 
associated spectroscopic strength [27]. 

Recently, experimental measurements of single-nucleon removal reactions including 7-ray detection have demon- 
strated that significant population of core fragment excited states may occur [27-34]. The inclusion of the core states 
within the theoretical framework [23,37], as foreseen in the original work of Sagawa and Yazaki [15], has lead to the 
proposal that such reactions (often termed "knockout" [35]) may be used as a spectroscopic tool [23,28]. To date, 
however, this approach has been largely confined to selected weakly bound halo and near dripline systems [28-33,36]. 

The objective of the present study was to undertake a systematic study of single-neutron removal reactions over 
a range of neutron-rich nuclei. For these purposes the light psd-shell nuclei were selected as the nuclear structure 
could be reliably calculated within the shell model. The region encompassed a number of nuclei of current interest 
and the production rates were relatively high. The use of high energy nucleon-removal reactions as a spectroscopic 
tool could thus be verified on a number of near stable nuclei with well known ground state structures. Additionally, 
the structural evolution with isospin across the psd-shell, as expressed in the core fragment observables, could be 
explored. As will be demonstrated, even inclusive measurements of the core fragments when executed using a broad 
range high acceptance spectrograph offer a means to survey changes in structure over a wide range of isospin in a 
single measurement. 

Earlier measurements of halo nuclei have suggested that the influence of nuclear and Coulomb dissociation on 
the core fragment longitudinal momentum distributions is relatively weak [12-14,38]. In order to explore further 
such effects, both Carbon and Tantalum targets have been used in the present work. The results obtained for the 
longitudinal momentum distributions and cross sections on the Carbon target have been briefly described elsewhere 
[39]. Here further details of the experimental techniques are given along with a detailed account of the theoretical 
models. In addition, the results obtained for the transverse momentum distributions using the Carbon target are 
presented, as are the longitudinal distributions from reactions on Tantalum. 

The paper is organized as follows. The experimental setup and techniques are described in Section II and the 
experimental results are presented in Section III. Sections IV to VII are devoted to the description of eikonal based 
modelling of nuclear and Coulomb dissociation. The results and comparison to calculations are discussed in Section 
VIII. A discussion, in the light of the present results, of the utility of single- nucleon removal as a spectroscopic tool 
is presented in section IX. The paper concludes (Section X) with a summary and perspectives. Explicit analytical 
formulae pertaining to the Coulomb dissociation calculations are presented in Appendix A. The results of the shell 
model and cross section calculations are tabulated in Appendix B. 

II. EXPERIMENTAL METHOD 

The experiment was performed at the GANIL coupled cyclotron facility. The secondary beams were produced via 
the fragmentation of an intense (<~ 1/iAe) 70 MeV/nucleon 40 Ar 17+ primary beam on a 490 mg/cm 2 thick carbon 
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target. The reaction products were collected and analysed (B/?=2.880 Tm) using the SISSI device [40] coupled to 
the beam analysis spectrometer. The resulting secondary beam was composed of 12 ~ 15 B, 14 ~ 18 C, 17 ~ 21 N, 19 ~ 23 
and 22_25 F nuclei with energies between 43 and 68 MeV/nucleon and intensities ranging from ^600 pps ( 15 C) to 
~1 pps ( 25 F). The energy spread in the beam, as defined by the acceptance of the beam analysis spectrometer, was 
AE/E—1%. Secondary reaction targets of carbon (170 mg/cm 2 ) and tantalum (190 mg/cm 2 ) were used. 

Owing to the large energy spread in the secondary beam, an energy-loss spectrometer was required to undertake 
a high resolution measurement of the core fragment momentum distributions [8]. In the present case, the SPEG 
spectrometer [41] was employed and operated at a central angle of zero degree in a dispersion matched mode for 
which an intrinsic resolution of 6p/p~ 4.5x 10~ 4 (FWHM) was achieved. The final resolution, including target effects 
was 6p/p=3.5xl0~ 3 (FWHM). The overall momentum acceptance of the spectrometer was 7%, which permitted the 
momentum distributions for the fragments resulting nuclei from one-neutron removal on all the nuclei of interest to 
be obtained in a single setting for each target (BpspEG—^-551 Tm for the carbon target and B ( ospbg=2.615 Tm for 
tantalum) . 

Importantly, the broad angular acceptances of the spectrometer — 4° in the horizontal (bending) and vertical planes 
— provided, in the case of the carbon target, for almost complete collection (see below) of the core fragments, obviating 
any ambiguities in the integrated cross sections and longitudinal momentum distributions that may arise from limited 
transverse momentum acceptances [42-44]. In the case of the tantalum target (Section VIIC and ref. [45]), the effects 
of multiple scattering and Coulomb deflection resulted in greatly reduced effective transverse momentum acceptances 
that curtailed the extraction of any reliable transverse momentum distributions or cross sections. An investigation of 
the effects of incomplete transverse acceptances on the longitudinal momentum distributions is presented elsewhere 
[45,46]. 

Ion identification at the focal plane of SPEG was achieved using the energy loss derived from a gas ionisation 
chamber and the time-of-flight between a thick plastic stopping detector and the cyclotron radio-frequency. Additional 
information was provided by the residual energy measurement furnished by the plastic detector and the time-of-flight 
with respect to a thin-foil microchannel plate detector located at the exit of the beam analysis spectrometer. Two 
large area drift chambers straddling the focal plane of SPEG were employed to determine the angles of entry of each 
ion and, consequently, allowed the focal plane position spectra to be reconstructed. The calibration in angle was 
performed using a tightly collimated beam and a calibrated mask placed at the entrance to the spectrometer. The 
momentum of each particle was derived from the reconstructed focal plane position. Calibration in momentum was 
achieved by removing the reaction target and stepping the mixed secondary beam of known rigidity along the focal 
plane. This procedure also facilitated a determination of the efficiency across the focal plane for the collection of the 
reaction products — the range of angles accepted by the spectrograph being restricted at the limits of the focal plane. 
Where necessary, corrections were then applied to account for any reduction in efficiency. 

The angles of incidence on target of the secondary beam particles were determined using two beam tracking detectors 
(each comprising four 10x10 cm 2 drift chambers) located in the analysis line of the spectrometer. The calibration in 
incident angle was derived from the trajectories reconstructed at the focal plane of the spectrometer for a measurement 
made with the target removed and the spectrometer set to the same rigidity as the analysis line. Consequently, the 
transverse momentum distributions for the core fragments could be reconstructed on an event-by-event basis from the 
incident projectile angle and the core fragment outgoing angle. The beam envelope was approximately Gaussian in 
form and characterized by half-angles of A^q/2^0.35 and A0!/ 2 =O.5 O (FWHM) in the horizontal and vertical planes 
respectively. Owing to the superior resolution in the determination of the angles in the bending (horizontal) plane 
of the spectrometer (68—0.1°, (50=0.4° (FWHM)), the transverse momentum distributions presented here have been 
reconstructed in this plane (p x ), with final resolutions including multiple scattering of some 5% being achieved. 

The intensities of the various components of the secondary beam were derived from a measurement of the primary 
beam current, which was recorded continuously during the experiment using a non-interceptive monitor, with respect 
to runs taken with the secondary reaction target removed and the spectrometer set to the same rigidity as the 
beamline. A redundant check was also provided by counting rates in two microchannel plate detectors placed in the 
incident secondary beam: one, as noted above, at the exit of the beam analysis spectrometer and another located 
on the upstream side of the secondary reaction target. The final cross sections were determined using an average of 
these three normalisations and the uncertainties quoted include contributions from both the statistical uncertainty 
and that arising from the determination of the secondary beam intensities (typically <~ 7%). 

The core fragment angular distributions were inspected to ensure that no events were lost due to the finite acceptance 
of the spectrometer (Af2=5 msr). For the reactions on the C target only in the cases of the broadest distributions 
did the losses exceed a few percent. The number of events not detected was estimated based on extrapolations of 
Gaussian adjustments to the measured angular distributions. Corrections were also applied where necessary to those 
nuclei falling near the limits of the focal plane (see above). The final uncertainty in the cross section includes an 
estimate (<~ 5-10%) of the uncertainties in these two corrections. As discussed in Section VIIC, for the reactions on 
the Ta target extremely broad angular distributions were encountered, which precluded any reliable estimates of the 
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transverse momentum distributions and cross sections to be made. 

In many instances — most notably 14 B, 15 ' 16 C and 17,18 N — asymmetric longitudinal momentum distributions 
exhibiting low momentum tails were observed (Fig. 1). The origin of these events is discussed in Section VIIIA. The 
cross sections reported here include these events (typically <~ 5%). 

To compare the measured distributions with the theoretical ones, all broadening effects inherent in the measurements 
should be taken into account. These effects arise from the differential energy losses of the projectile and the fragment 
in the target, energy and angular straggling in the target and the detector and spectrometer resolutions. In addition, 
for the longitudinal momentum distributions the Lorentz transformation from the laboratory to projectile frame of 
reference must be taken into account. In order to provide an estimate of the relative importance of these effects on 
momentum distributions, a Monte Carlo based simulation was developed. As an example, an evaluation of the effects 
for two nuclei with distributions representative of those encountered here are provided in Table I. The overall effect 
in the case of the longitudinal distributions is some 10-13% and is dominated by the Lorentz contraction. In the case 
of the transverse momentum distributions, the broadening is relatively weak (at most some 5%) and is dominated by 
the angular resolution of the spectrometer. 

The widths of the momentum distributions were derived from Gaussian adjustments to the central region (~FWHM) 
of each distribution, thus avoiding any bias introduced by low momentum tails. The use of other lineshapes (such as 
a Lorentzian), or a simple statistical analysis [14] produced essentially identical results. The widths, noted FWHM cm 
in Table II, are quoted in the projectile frame and have been corrected for the various broadening effects discussed 
above. 



III. RESULTS 

The core fragment longitudinal momentum distributions measured using a carbon target are displayed in Fig. 1 . In 
order to facilitate their comparison, each distribution is displayed over the same total momentum range of 700 MeV/c. 
The estimated widths, taking into account the various experimental effects discussed above, are shown in Fig. 2 and 
are listed in Table II (FWHM™). The corresponding single- neutron removal cross sections are tabulated in Table II 
and the evolution along the isotopic chains is presented in Fig. 3. 

The transverse momentum distributions p x from breakup on the carbon target are presented in Fig. 4. The same 
total momentum range (here 600 MeV/c) has been used to display all the results to facilitate the comparison. The 
widths extracted from the measured distributions, taking into account all broadening effects (Section II), are listed 
in Table II (FWHM™). A comparison of the widths of the longitudinal and transverse momentum distributions 
is provided in Fig 2. The transverse distributions are systematically somewhat broader, a feature already observed 
in reactions with stable beams [47]. More interestingly, the transverse distributions exhibit the same trends as the 
longitudinal distributions, suggesting that the sensitivity to the projectile structure expected for the latter is also 
present to a similar degree in the transverse distributions. As will be discussed further in Sections V-B and VIII, 
the transverse momentum distributions present somewhat more complex forms than the longitudinal distributions. It 
should be stressed that the widths quoted here are only meant to serve as a comparative guide. 

For a number of nuclei ( 14 B, 15 ~ 17 C and 17-19 N) the core fragment longitudinal momentum distributions from 
reactions on a Ta target were also measured (Fig. 5). The momentum widths, which are seen to be almost identical 
with those obtained on the C target, are listed in Table III. As outlined in Section II, the corresponding transverse 
momentum distributions were observed to be much broader than the acceptances of the spectrograph (see, for example, 
Fig. 16). As such no reliable cross sections could be derived. 

As noted in Tables II and III, a number of the nuclei studied in the present work have been investigated elsewhere. 
In the cases of 14 B and 15 . 17 > 18 C measurements made at similar energies on a Be target by Bazin et al. [48,49] found 
longitudinal momentum distributions with widths in good agreement with those presented here. The associated 
single-neutron removal cross sections are, however, significantly smaller than those measured in the present work. As 
pointed out in our earlier paper [39] the origin of this discrepancy lies in the rather limited acceptances of the A1200 
fragment separator. This is clearly apparent from the transverse momentum distributions presented here (Fig. 4) [45] 
and is also confirmed by more recent measurements of 14 B [30] and 16 ' 17 C [32] undertaken using the high acceptance 
S800 spectrograph (Table II). 

Very recently a measurement for 16 C (83 MeV/nucleon) on a C target has been reported [50]. The results, which 
were obtained using a new time-of- flight technique to deduce the core momentum distribution [51], are in very good 
agreement with those reported here. 

A measurement has also been carried out with the timc-of-flight technique, at an energy somewhat higher 
(72 MeV/nucleon) than that employed here, of the breakup of 23 by carbon [51]. The resulting 22 longitudi- 
nal momentum distribution is in good agreement with the present work. Moreover, as discussed in Section VIII-D, 
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the corresponding one-neutron removal cross section (233±37 mb) is in good agreement with that calculated here 
(Tables II and IV) 1 . 

At much higher energies (^900 McV/nucleon), the single- neutron breakup of 17 C by carbon has also been measured 
[38,52]. A width (FWHM) of 141±6 MeV/c, slightly broader than reported here, was extracted [38] together with 
a cross section of 129±22 mb [52]. Recently an experiment using the same setup as that of refs. [38,52] has been 
carried out to explore single- neutron removal on neutron-rich isotopes of N, O and F at some 900 MeV/nucleon. The 
preliminary results for the momentum distributions — in particular for 23 O (Section VIII-D) — are in very good 
agreement with those reported here [53]. 

In the case of reactions on the Tantalum target (Table III), there is a relative paucity of work with high-Z targets 
with which comparison can be made. Indeed, of the nuclei measured here, the literature reports results for only 14 B 
and 15,17 C. In the case of the former, very good agreement is found with the measurement made employing the S800 
spectrometer [30] , whilst the earlier work of Bazin et al. found a somewhat narrower core momentum distribution [49] . 
The core momentum distribution measured for 15 C in the same experiment [49] is in reasonable agreement with that 
observed here. In the case of the heavier isotope, 17 C, a preliminary study [80] reported a momentum distribution 
consistent with that found in the present work. Unfortunately the more complete study of Maddalena and Shyam 
[79] does not quote any widths for the momentum distributions. 

In order to examine quantitatively the relationship between the projectile structure and the measured distributions 
we now turn to a detailed development of the necessary reaction theory. 

IV. EXTENDED GLAUBER MODEL FOR NUCLEAR DISSOCIATION 

The ensemble of data presented in the previous section constitutes a test not only for the structure models of 
neutron rich nuclei in p-sd shell but also for the description of the reaction mechanisms involved. The aim of this 
section is to provide a model that, starting from realistic projectile wave functions and taking properly into account 
the reaction mechanism, can explain the momentum distributions and cross sections. In the formal development, the 
principle features of which are similar to approaches developed by Esbensen [18], Bertsch et al. [20,54], Negoita et al. 
[22] and Tostevin [23,37], we have attempted to retain the role played by the wave function via its Wigner transform 
as this reveals clearly the momentum content. The formulae obtained for the fragment momentum distributions and 
cross sections explicitly display, in the spirit of the Glauber model, the distorting functions arising from the reaction 
mechanism. The details concerning the practical calculation of some basic ingredients — in particular the S-matrix 
elements — are also provided. 

We shall assume that the ground state of the projectile ( J T ) could be approximated by a superposition of configura- 
tions of the form [I£ £g) nlj] J , where I£ denotes the core states and nlj are the quantum numbers specifying the single 
particle wave function of the last neutron, taken here as Woods-Saxon wave functions evaluated using the effective 
separation energy S%?* = S n + E^ x (E ex being the excitation energy of the core state) . We neglect coupling of core 
states to the final state and dynamical excitation of excited core states in the reaction. In this approximation the 
reaction can populate a given core state only to the extent that there is a nonzero spectroscopic factor C 2 S{I^, nlj) in 
the projectile ground state. When more than one configuration contributes to a given core state, then the total cross 
section for one neutron removal is written, following refs. [23,37], as an incoherent superposition of single particle 
cross sections: 

v-mil*) = £ C 2 S(i:, nlj) a sp (nlj, S%"). (1) 

nlj 

The total inclusive one-neutron removal cross section (cr?^ tt6er ) is then the sum over the cross sections to all core 
states. A similar relation holds for the momentum distributions. The coordinate system used in the calculations is 
sketched in Fig. 6, whereby the impact parameters for the neutron and the core are given by, 

K = R ± + (2) 

I = R ± - -^s. (3) 



x As no reliable beam intensity normalization was available in the present study for 23 0, no experimental cross section could 
be extracted 
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We shall neglect recoil effects, so that (A c 3> 1), 

b n = b + s + o(^j , (4) 

tc = *+o(±-). (5) 

Since the core states are not coupled by the interaction, the core plays a spectator role. Thus it is sufficient to 
consider only the neutron degrees of freedom as described by the wave function ^jm corresponding to the coupling 

scheme [[I <g> i]j <8> I C ]J ■ We assume only one bound state in the projectile and we need to consider density matrix 
elements of the form: 

PM'm(^) =< *JM'|*JM >spins, (6) 

where <> sp » ns means average over spin coordinates. A little angular momentum algebra leads to, 

P J M' M (r) = E(-) J " 2j " /c " l ^ C oooGM^W (jljl, ±a) W (JjJj, I C X) R 2 nl3 {r)Y x ,{f), (7) 
with the property, 

Tr P = "i E Pm>m 5 m>m = R 2 nlj {r)/A-K = Pva i(r), (8) 

J M 

with R n ij(r) the radial part of the single particle wave function. It is useful to introduce also the projected density, 

/oo 
dz PM'M(r), (9) 
-OO 

Pvai(s) = i Yl Pmm{ s )- (!0) 

J M 

If S c and S n are the S-matrices in impact parameter representation for the core and neutron-target interactions, 
the absorption cross section (or stripping) is given by, 

(Tabs = i / ^SSl < <t>k\^abs\^JM > | 2 , (11) 
•* J k M 

where 4>k are scattering states and fi 2 bs = (1 — |S n | 2 )|S c | 2 is the transition operator for neutron absorption. As the 
neutron is absorbed, only the scattering states are available to the core. In this case the closure relation is, 

&|=1. (12) 

fe 

Combining Eqs. (11) and (12) we obtain, 

^abs = Y 2 £ I s j df i 1 ~ \ S ^+ l^(&)| 2 flU(0> (13) 

which we rewrite in the form 

(Jabs = J dbds{\ - \S n (b + s)\ 2 )\S c (b)\ 2 p val (s) = J dsD abs (s)~Pvai(s), (14) 

where D abs is the distortion kernel [22] . This kernel has a very intuitive physical interpretation: it is the convolution 
product of the survival probability for the core and the probability for the absorption of the neutron, as originally 
defined by Hiifner and Nemes [6]. In fact the general absorption operator can be decomposed as follows, 

1 - |S C | 2 |S„| 2 - |S C | 2 (1 - \S n \ 2 ) + \S n \ 2 (l - \S C \ 2 ) + (1 - |5„| 2 )(1 - |5 C | 2 ). (15) 
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where the first term corresponds to neutron absorption, the second to core absorption and the double scattering term 
to the absorption of both the core and the neutron. For the last two (inelastic) channels the total cross section is 
formally identical to Eq. (14) with appropriate redefinitions of the distorting kernels. 

For diffraction a similar formula to Eq. (11) holds, except that the transition operator is replaced in this case by 
Qdiff = S c S n — 1. We shall again assume a structureless continuum and treat it via sum rules. In addition it is 
assumed that the projectile has only one bound state. Since the scattering states should be orthogonal to the ground 
state, the closure relation is in this case: 



\^JM >< + E l^fc ><M = 1- 



(16) 



M 



Clearly, the "one bound state assumption" will lead to an overestimation of the breakup cross section, since any 
additional bound state will subtract cross section [87]. The orthogonality condition allows us to replace the "—1" in 
the definition of the diffraction transition operator by any function which does not couple the nucleon coordinates. 
For the purpose of convergence, the most convenient form for this operator is Cldiff = S c (S n — 1). With this, the 
total cross section for diffraction is, 



' ff = ~h / * E < ^JM\\^diff\ 2 \^JM > ~ I < ^JM\£ldiff\^JM> >| 2 , (17) 

J •* Km m'm J 



which demonstrates that the diffraction cross section is given by the fluctuation of the transition operator in the 
ground state. In terms of the density matrix (6-10) the nuclear diffraction cross section is, 



Vdiff 



J2 



dbds \S c (b)\ 2 \l-S n (b + s)\ 2 p val (s) 

2 



M'M 



db 



dstt d i ff {b,s)p J M , M {s) 



(18) 



The first term is similar in structure to (14) and provides the main contribution to the diffraction cross section. The 
second is a small correction arising from the orthogonality requirement. One can further simplify this correction by 
observing that it arises essentially from the diagonal elements of the density matrix. Indeed, if one considers the 
following integral, 



J dsdz fi (b, s) p J M , M (r), 



and expands the density in multipolcs, 



P J M'M(r f )=J2px„(r)Px„(8)e i ' 



(19) 



(20) 



where cos# = z/r and <p — 4>{ x , y) an d assume that the dependence of O on angles of s is weak and may be neglected, 
we are left with, 



(21) 



Therefore the main contribution to such integrals comes from multipoles with p = and only the diagonal elements 
of density matrix contribute (M' = M, see (7)). Our final formula for diffraction cross section is then, 



a diff = J dbds\S c (b)\ 2 \l - S n (b + s)\ 2 p vai (s 
/ d b\Sc(b)\ 2 I ds(l-S n (b + s))p J MM(s) 

M J 

As mentioned earlier the first term is dominant and can be written in the form, 



CT *// = j d * D diff(s) Pval(s), 



(22) 



(23) 
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with D diff ~ D abs . 

Further simplifications arise if one observes that for spin independent transition operators one can neglect the 
intrinsic spin of the nucleon and the coupling to core spin. In this case the total spin J is replaced by the angular 
momentum I. Thus, for example the density matrix (7) becomes, 



I 2 



(24) 



which has the same properties as in Eqs. (8-10). The total cross sections for absorption and diffraction become, 

a abs = J dbds (l - \S n (b + s}\ 2 ) \S c (b)\ 2 p val (s), (25) 



Vdiff 



= Jdbds\S c {b)\ 2 \l-S n (b + s)\ 2 p val {s) 

-^E/ d b\Sc(b)\ 2 J ds (l-S n (b + gj) p l mm (s) 



(26) 



Finally one may note that the absorption cross section is fully equivalent with the corresponding cross section 
obtained by Hussein and McVoy in the core spectator model [19]. 



V. MOMENTUM DISTRIBUTIONS 



A. Longitudinal momentum 

As final state interactions are neglected the scattering states are taken as plane waves. The basic matrix element 
which gives the localization probability in momentum space is, 



^ = tAtT2 E \{e- rk % m M)XiMAm(b,mJM(r,a,0) 2 
dk (2tt)3J2 MMims I 



(27) 



Integration is taken over nucleon (r, a) and core spin (£) coordinates. The intrinsic momentum distribution (Wo(k)) is 
obtained by choosing fi = 1. For absorption and diffraction one uses the appropriate operators defined in the previous 
section. After applying some angular momentum algebra, one finds, 



dP 



dk (2tt) 3 



12 

1 m, 



dre^nib^R^Y^if) 



(28) 



The longitudinal momentum distribution (fe z ) is obtained by integration over the unobserved components (k x and 
k y ). To obtain closed formulae it is useful to introduce the partial Wigner transform of the wave function, 



w m (s,k z ) = —L= [ dze tk * z Ri{r)Y lm (r), 
\Z2ttP J 



in terms of which the total Wigner transform is, 

W(s, k z ) = ^2 \ w m{s, k z )\ 2 , 



(29) 



(30) 



with the properties, 



J dk z W(s,k z ) = p va i(s), 
J dsdk z W(s,k z ) = 1, 
J dsW(s,k z ) = W (k z ). 



(31) 
(32) 
(33) 
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The longitudinal momentum distribution for absorption is thus calculated as, 

( Jr) = / Sdg i 1 ~ ]Sn{i+ • <r)|2 ) |5c(fe)|2 w{g > kz) = J dgD -bs(s)w(s, k z ). 

For diffraction the situation is somewhat more complicated and three terms must be considered, 



da N 


f 


dk Z/ 


' diff 


da N 


r 


dk Z/ 


diff 


da N 


r 


dk Z/ 


diff 



dbds\S c (b)\ 2 \S n (b + s}- 1| 2 W(s,k z ), 
£ db\S c (b)\ 2 ds w mi (s,k z )w^, l2 (s,k z ) 

mirii'm-2 

x f ds! p mim (si) (s n (6 + si)-l) j ds 2 p* m2m {s 2 ) (s*(b + s 2 )-l) 
-23? ]T f db\S c (b)\ 2 f dSi p mim (S*(b+ Si) - l) 

mmi J J 

x j ds 2 w m {s2,k z )w* mi {s 2 ,k z ) (s n {b + s 2 ) - l) . 



(34) 



(35) 



(36) 



(37) 



To a good approximation one can use again the fact that the main contribution to diffraction arises from the diagonal 
part of the density matrix. In this case the second and third terms become, 



da 

dk z 

da 

dk z 



(2) 

diff 
(3) 

diff 



I db\S c (b)\ 2 f ds\w m {s,k z )\ 2 j d3i(s n $+3{)-l)p mm {8{) 

m 

-23? ]T [db\S c (b)\ 2 f dSi/5 rom («i)(<S:(6 + ?i)-l) 
x J ds 2 \w m {s 2l k z )\ 2 (s n {b + s 2 ) - 



(38) 

It can be checked that integration over k z , in the range (— oo,+oo) leads exactly to the integrated cross sections 
derived in the previous section. It has been shown by Bonaccorso and Bertsch [55] that the error introduced by these 
limits is less that 5 % for the beam energies considered in the present experiment. 



B. Transverse momentum 

Transverse momentum distributions (k x , k y ) are obtained by projecting the probability (28) onto the axis of interest. 
For example we have, 



dP 

dk x 2ttI 



— d V dz / dxe ik " x 



Sl{b,s)Ri{r)Y lm {r) 



(39) 



For the other transverse direction one simply exchange the indices x < — ► y. The corresponding cross sections are 
given by 



da 

dk. 



x / abs 



7^pY,J dbn[l-\S n (b n )\ 2 ] J dydz J dxe ik * x S c (\b n - 3\)Rt(r)Y lm { 



(f) 



and 



£-) = At E / d b\Sc(b)\ 2 [ dydz f dxe ik * x (s n (\b + Xs\) - l) Ri(r)Y lm (r) 



(40) 



(41) 



with A = A as a measure of recoil effects. In the derivation, we assume small excitation energies in a structureless 
continuum (plane wave approximation) and neglect the orthogonality requirement described in the preceding subsec- 
tion. If A = 1, integration of equations (40-41) over k x leads to the total cross sections (14) and (23) respectively. It 
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should be noted that the transverse momentum distributions are symmetric around k x = 0. There is no asymmetry 
in the nuclear breakup transverse momentum distributions, principally due to the straight line approximation for the 
trajectory and the inherent neglect of the conservation of the energy. This symmetry means that only half of the 
distributions need to be computed. The necessity to evaluate a five dimensional integral for each value of k x , for two 
different transition operators and for each core state renders the calculations onerous. One may reduce the complexity 
by assuming that the angle dependence (b, s) of the transition operators is weak and can be replaced by an average 
value. For example, in the case of diffraction one has, 

S n (b, s) = ^ J d(j) [S n + A 2 s 2 + 2Xbs cos </>) - l) . (42) 

The use of average transition operators allows a straightforward integration over angles of b and reduces the integration 
domain to (0,oo) for x,y and z variables. The computational time is thus reduced by one order of magnitude and 
the precision in the calculation is significantly increased. We have checked that averaged transition operators contain 
almost the same transverse momentum components as the original ones, at least for impact parameters in the range 
of the strong absorption radius. 

Equations (40) and (41) illustrate the essential difference between the longitudinal and transverse momentum dis- 
tributions. The former arc given essentially by the Wigner transform of the valence nuclcon wave function weighted 
by distortion kernels for absorption and diffraction, whilst for transverse momentum distributions additional compo- 
nents appear due to nuclear interaction. The neutron-target interaction influences strongly the transverse momentum 
distribution for diffraction and leads to a broader distribution as compared to that arising from absorption (see 
Fig. 14). 

VI. COULOMB DISSOCIATION 

In this section we describe the Coulomb dissociation within the framework of the eikonal approximation which is 
formally equivalent to a first order perturbation theory. We use the long wavelength approximation for the transition 
operator and obtain a general formula for any electric multipolarity {Elm). However, the calculations are done only 
for dipole and quadrupole transitions, since these give rise to the principle contributions to the cross sections. 

The Coulomb excitation amplitude for a projectile in the field of a target may be expressed in terms of the electric 
multipole matrix elements characterizing the electromagnetic decay of nuclear states. If the charge distributions of 
the two nuclei do not overlap during the collision, then the relative motion takes place along a classical Rutherford 
trajectory. At very high energies, the trajectory is well approximated by a straight line and the first order eikonal 
approximation should give reasonable results for momentum distributions and integrated cross sections. 

Our starting point for the Coulomb amplitude (fc), is the result of Bertulani and Baur [56] in a first order eikonal 
approximation (formally equivalent to the first order Born approximation), 

fc(Q,q) =i^R 2 Y,i m i(^) 1 G Elm <P m (Q) M(Elm), (43) 

Ira 

where M(Elm) are the matrix elements for electric transition of multipolarity (Im) 

M(Elm) = Z 3 f Y lm(fj) <PiW df, (44) 

and ¥>»(/) are initial (final) state of the projectile. This is the long wavelength limit of the transition operator written 
in terms of cluster coordinates fj(j = 1,2). The two clusters are characterized by masses (m,j), charges (Zj), and 
momenta (kj). The relative motion of the outgoing clusters is described in terms of the relative momentum, 

q = (m 2 ki - mife)/(mi + m 2 ), 

whilst the momentum change in the scattering (Q) is given by Q — 2k sin | , where 9 and k are respectively the 
scattering angle and the incident momentum in the center of mass. The other notations in formula (43) are: the 
target charge, Z t ; the velocity of the projectile in units of the speed of light, /3; the relativistic Lorentz factor, 7; the 
fine-structure constant, a; and the interaction radius, R. In addition 

rko = e + E q = e + — — , 
2fj, 
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is the excitation energy, i.e. the sum of the (absolute) binding energy and the kinetic energy of the separated clusters. 
The relativistic functions of Winthcr and Adler [57] are used in the form Gei™ = i l+m G(f3,j) with the complex 
phase factorized out, whilst the nondimensional functions $ m (Q), also defined in refs. [56,57], contain information on 
reaction mechanism whereby 



J m {QRx) K m ( \ xdx, 



and standard notation has been employed for the Bessel functions. The transition matrix elements (44) are calculated 
using a simple shell model wave function for the ground state, 

<Pi(i f )=K(r)Y LM (f). (45) 

Note the change in notation with respect to the previous sections for the quantum numbers of the wave function: 
L, M instead of I, m. If final state interactions are also neglected, the continuum may be considered to be structureless 
and the final state wave function may be treated as plane waves. This is a good approximation for small excitation 
energies. 

Under the above approximations, the matrix element in (44) is given by, 



M(Elm) = (Zrfi - {-) l Z 2 fy j dre^r 1 Y lm {r) K L (r) Y LM {f) 
= V^zf f It ^r^A" 1 ^) C^L / r*drr l j x ( q r)1l L (r), (46) 

with an obvious notation for the effective charge Z e ^ and /?i(2) = ^2(1) /{ m i + 1^2}- Since the spin orientations are 
not specified, the differential cross section for Coulomb excitation is obtained by averaging the square of the Coulomb 
amplitude over the magnetic projections, 



M 



The main contribution in (46) is given by dipole (El) and quadrupole (E2) transitions, therefore only terms with 
I = 1 or I = 2 are included. The complexity of Eq. (47) arise from the evaluation of integrals, 

/>oo 

/m ima (0 = 2R 2 e / QdQ$* mi (Q,0$ m2 (Q,0, (48) 
Jo 

where, 

_ ujR _ uj R 
jv c /fy 

is the adiabaticity parameter, defined in terms of the excitation energy (%uS) and the minimum impact parameter R, 
given by 

R = R p + Rt + 



±E labl ' 

which includes a correction due to the deviation of the trajectory from a straight line [57] . The function (48) is obtained 
in general by numerical integration except for the diagonal term which admits a simple analytical expression. However, 
one can profitably perform first the integration over the azimuthal angle (</>) of the relative momentum, 

dq = q 2 dq sin 9d9d(f>. 

This integration automatically selects only the diagonal terms in (48) and, 

fmm — Km{£)i 

with functions K given by: 

K,(0 = (Kl-Kl), 
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K 1 {£,) = -^K Q K 1 -(Kl-Kl), 



K 2 ^) = ^Kl + {Kl-Kl), 



where K m are modified Bcssel functions of the first kind and K _ m = K m . After some simple calculation we obtain 
the following expression for the differential cross section: 



// 7 eff 



q 2 dqsin9d6 ^/Att ^ ' vc/ r 

x Zf Zf Gi im Gi 2Tn K m (£) 



x ihhxAq) iu 2 \ 2 (q) Cqoo 1 C^q 2A 2 C^q 2 q s C^i^ jf^Ai} Y so(q), 

where Iu\(q) is a shorthand notation for the radial integrals appearing in (46). The summation runs over all 
quantum numbers except L. The angular dependence of the cross section is given by functions Y$o which are the 
standard spherical harmonics defined without the phase factor e tm ^ (where <f> has already been integrated). The 
energy dependence of the reaction mechanism is governed by the functions G and K. However, the magnitude of the 
cross section for a given multipolarity transition is mainly determined by the effective charge. If the clusters have 
equal charge to mass ratio, then the dipole transition has vanishing cross section in this approximation. This is readily 
understood from classical arguments since in this case the dipole field acts on the two clusters with the same force 
in the same direction and does not lead to dissociation. This is a consequence of the assumption of a well defined 
cluster structure for the projectile. Experimentally an appreciable, but not complete suppression of the El transition 
occurs. The interference term E\E2 does not contribute to the total cross section. 

Various observables may be readily obtained from (49) by appropriate integration or change of variables. The total 
Coulomb dissociation cross section ac = &ei + &E2 can be obtained by numerical integration of equation (49). In 
practice we have used the closed formulae given in Appendix A. For momentum distributions a change of variables is 
made such that, 

dq = dq x dq y dq z = 2nq 2 dq sin 9 d9 = 2tt q r dq r dq z . 
The radial momentum distribution is obtained by integrating over the longitudinal momenta, 

r i n r ta , (<l 



d< , £ *• E W» (f ) - 2 E 1°° *Mr» v , 

a, b a, b — even 



dq, 



where we have used a shorthand notation for the general cross section (49). In addition, q 2 = q 2 + q 2 , cos 6* = q z /q, 
sin# = q r /q and (a) denotes all summation indices appearing in (49), different from S. There is no E1E2 asymmetry 
in the radial momentum distribution since the interference term in (49) contains only odd S-values. Similarly, one 
can demonstrate that the longitudinal momentum distribution takes the form, 



da 

dq z 



(q z > 0) = J2 r A S^ Y ^ (-) « d 1> 
a ,S J 4* V q ' 



and 



This leads imediately to, 



-^(fo < 0) = £ / V) S A^q)Y S o ( q ~) q dq. 



—J— — \1z < 0) = — (q z > , 

dq z dq z 
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d(JEiE2, ^ nN da E \Ei, ^ n , 

— (q z < 0) = q z > . 

dq z dq z 

From the above relations it can be seen that the loss of information — concerning e.g. the E1E2 interference term 
- inherent when an integration over all variables is performed, may be partially compensated for by measuring 
the longitudinal momentum distribution. In an inclusive measurement, such as that performed here, only the core- 
like particle momentum is measured. To compare with the data, we have to transform the theoretical momentum 
distribution which is given as a function of relative momentum to a function of fragment momentum. This is done 
most easily in the projectile rest frame, taking into account momentum conservation, 

q = Mi - A;92 
<Zi + q\ = 0, 



and 

da 



da_ f 

dq z i J 



dq z dq z2 f(q z ) S{q z - j3\q z i + (3 2 q Z 2) S(q z i + q z z) = f(q z i)- 



where / is a generic notation for the theoretical momentum distribution. A similar formula holds for the other 
fragment. 

VII. S-MATRICES AND OPTICAL MODEL POTENTIALS 

The remaining physics is to describe interaction of the core and the removed nucleon with the target. These enter 
through the associated S-matrices, S c , n , expressed as a function of impact parameter. Previously, S-matrix calculations 
have been based on the optical limit of the cikonal model [20,22,23,28]. In this approach the nucleus- nucleus phase 
shifts arc entirely determined by nucleon-nucleon collisions in the density overlap volume. The energy dependence is 
dictated by the total nucleon-nucleon cross sections vnn- However, diffraction dissociation is sensitive to the refractive 
power of the optical potential and these effects are not very well controlled in the above approximation. 

Recently, a more fundamental approach has been introduced by Bonaccorso and Carstoiu [58] and Tostevin [37], 
whereby the G- matrix interaction of Jeukenne, Lejeune and Mahaux (JLM) [59], which is obtained in a Brueckncr- 
Hartree-Fock approximation from the Reid soft core nucleon-nucleon potential, has been employed. This interaction is 
complex, density and energy dependent and, therefore, provides simultaneously both the real and imaginary parts of 
the optical potential. The optical potential calculation with this interaction is described in detail in refs. [58,60]. The 
single particle densities for the core and target were generated in a spherical Hartree-Fock + BCS calculation using 
the density functional of Beiner and Lombard [61]. The strength of the surface term in the functional was slightly 
adjusted in order to reproduce the experimental binding energy for each nucleus. The matter radii resulting from 
these calculations are in good agreement with the existing experimental values of Liatard et al. [62] and with results 
of relativistic mean field calculations (RMF) by Ren et al. [63-65], as shown in Fig. 7. In Fig. 7, other experimental 
radii values derived from reaction cross sections at higher energy are also presented [66,67]. Except for the case of 
oxygen isotopes, these values are systematically slightly smaller than those obtained by Liatard et al.. 

The resulting optical potentials were renormalised in order to reproduce the total cross section for neutron target 
interaction in an eikonal calculation including noneikonal corrections up to second order. For the core target potentials 
the renormalisation constants have been taken from ref. [60]. It should be noted that the potentials are strong and 
the eikonal expansion (as defined by Wallace [68]) does not converge at low energies. 

The resulting S-matrix elements and transmission coefficients are displayed in Fig. 8 for n+ 12 C and 14 C+ 12 C at 
energies of 30 MeV and 30 MeV/nucleon respectively. One sees clearly substantial changes in the shape and magnitude 
of all matrix elements for the neutron if higher order noneikonal corrections are taken into account. These effects are 
most pronounced at low impact parameters where neutron absorption profile shows an important contribution from 
trajectories reflected inside the barrier superimposed on a characteristic strong absorption at the nuclear surface. The 
eikonal approximation in lowest order underestimates the interaction range and the absorption in the nuclear interior. 
Given the surface dominance of the breakup reactions, this approximation will lead to an underestimation of the 
stripping and dissociation cross sections. The distortion kernels obtained from these S-matrix elements are displayed 
in Fig. 9. It is clear that diffraction dissociation is most sensitive to the noneikonal corrections. The asymptotic 
behaviour of the distorting kernels is most affected and this has important consequences for the calculation of the 
momentum distributions, since large impact parameters probe the low momentum content of the projectile wave 
function. 
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Finally, the neutron S-matrix has been checked against known experimental total cross sections for four targets 
[69]. The results are displayed in Fig. 10. The comparison is reasonably good for all targets except tantalum where 
the absorptive potential is too strong to be treated in the eikonal approximation. Nevertheless, we have been able to 
find normalization constants which reproduce, at least qualitatively the experimental cross sections. The second order 
eikonal calculation for n+ 12 C match reasonably well recently evaluated data for elastic and reaction cross sections in 
the range 20-100 MeV [70] (Fig. 10b). 

VIII. DISCUSSION 

A number of features are apparent from the systematics of the core fragment momentum distributions and associated 
single- neutron removal cross sections for reactions on Carbon (Figs. 1-3). Firstly, the crossing of the N=8 shell 
and N=14 subshell closures are associated with a significant reduction in the widths of both the longitudinal and 
transverse momentum distributions for 14 ' 15 B, 15 C, 23 and 24 > 25 F compared to the neighbouring less exotic nuclei. 
Such behaviour is a clear indication of the role played by the structure of the projectile owing to the large valence 
^2si/2 admixtures expected in the ground state wave functions (as discussed in detail below). Secondly, as a result 
of this and the weak binding of the valence neutron 14 B (S„=0.97 MeV) and 15 C (S„=1.22 MeV) exhibit enhanced 
cross sections in comparison to the neighbouring isotopes, suggesting a spatial delocalisation of the valence neutron 
orbital [39]. 

In order to analyse quantitatively the measurements presented here, we now proceed to make a comparison with 
the results of calculations using the extended Glauber model described above coupled with the results of shell model 
calculations. 



A. Cross sections and longitudinal momentum distributions: Carbon target 

The spectroscopic factors C 2 S{I* ,nlj) employed here were calculated with the shell-model code OXBASH [71] 
using the WBP [72] interaction within the lp-2sld model space. In order to avoid energy-shift effects [73] only pure 0, 
1 or 2 hui excitations were considered. The older Millencr-Kurath [74] interaction (PSDMK) was also investigated and 
the results were found, for the nuclei examined here, to be comparable to those obtained using the WBP interaction. 
Where known, the experimentally established spin-parity (J 7r ) assignments and core excitation energies [75-78] have 
been used for calculation of cross sections and momentum distributions. In all other cases, the shell-model predictions 
were employed. It was found that the cross sections were relatively insensitive to the excitation energies of the core 
states. A detailed listing of the spectroscopic factors and calculated cross sections are given in Appendix B 2 . 

Aside from the spectroscopic factors the principle parameters entering into the calculations were the S-matrices 
and the geometry of the Woods-Saxon potential used to define the single-particle wave function. As described in 
section VII, the JLM calculations of the S-matrices were verified through comparison with measured total, elastic and 
reaction cross section data. The strength parameter of the single-particle potential was fixed by fitting the known 
experimental one-neutron separation energy. The radius and diffusivity of the Woods-Saxon potential were fixed at 
r„, s =1.15 fm and 0^,8=0.5 fm for the isotopes of B and C, and r„, s =1.2 fm, and a ws =Q.d> fm for N, O and F. These 
values were chosen to provide a good global agreement with the measured cross sections. Better agreement is obtained 
if for example radius and diffusivity parameters are tuned locally. The sensitivity to the choice of these parameters 
may be illustrated by the example of 16 C breakup at 50 MeV/nucleon, for which a change in geometry to r„, s =1.20 
fm and 0^=0.65 fm leads to an increase of about 27 % in the total one-neutron removal cross section. As might be 
expected the shape and width of the momentum distributions were found to be rather insensitive to the radius and 
diffusivity since these parameters mainly affect the single-particle asymptotic normalization coefficient of the wave 
function. It should also be noted that in the present calculations the ground and excited states of the core were 
assumed to have the same density distributions. As such the same Woods-Saxon geometry was employed for all core 
states of each projectile. 

In order to facilitate the comparison of the calculated and measured momentum distributions, the former were 
filtered through a Monte Carlo simulation (Section II) to take into account the experimental broadening effects. As 
may be seen in Figs. 1 and 3, the measured distributions and cross sections for all the nuclei included in the present 



2 Since our original publication [39,45] the 22 ' 24 F and 23 calculations have been revised. 
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study, including those nuclei with well established structure, are well reproduced with the exception of 22 F, where the 
cross section is somewhat underestimated. 

It is, however, apparent that for a number of nuclei the calculated momentum distributions are slightly broader 
than the experimental ones 3 , in particular the low and high momentum wings are somewhat more pronounced. This 
effect appears to arise from the specific shapes of the distorting functions (-D(s), Section V) at low impact parameters, 
as shown in Fig. 9. The use of distortion kernels calculated with less realistic black disk S-matrix elements leads to a 
strong suppression of the high momentum components in the wave function, and momentum distributions consequently 
become much narrower. 

To illustrate the contributions from the different mechanisms leading to the removal of the neutron, the various 
contributions — absorption, diffraction and Coulomb dissociation — to the breakup of 14 B, 15 C, 17 C and 21 are 
displayed in Fig. 1 1 . The Coulomb dissociation cross section is typically less than 1 mb in all but the most favourable 
cases — 14 B and 15 C — for which the Coulomb induced breakup was estimated to amount to some 7 mb. As 
expected [17], absorption and diffraction result in distributions with very similar lineshapes. The diffraction cross 
section is, however, smaller than that for absorption for well bound states, whilst the two are essentially identical for 
weakly bound states (Fig. 3). This evolution with binding is also illustrated in Fig. 12 for single-neutron removal 
at 50 MeV/nucleon from an A = 17 system comprising a core and a single- valence neutron. In the case where the 
neutron occupies an s-wave configuration (left panel) the absorption and diffraction cross sections are almost equal 
independent of the binding energy. If the beam energy is increased, absorption becomes the dominant process, whilst 
at lower energies diffraction is favoured. For a d-wave valence neutron the cross section is dominated for all energies 
by absorption and the contribution from diffraction decreases as S n increases. 

The effects of applying the noneikonal corrections described in Section VII are also displayed in Fig. 12. These 
corrections lead to an increase in the total (absorption and diffraction) cross section: for example, for a valence neutron 
binding energy of 1 MeV, an increase of some 12% occurs. As expected [91] the effect is even more pronounced at a 
lower energies (e.g., some 19% at 30 MeV/nucleon). 

Following Brown et al. [89] we introduce a quenching factor R s — <J-i n /o'-i% lber in analogy with the z factor of 
Pandharipande et al. [90]. Individual values are plotted in Fig. 13 as a function of projectile mass number. Averaging 
over the 22 one-neutron removal reactions for which cross sections were measured here, one obtains R s = 0.98 ± 0.16. 
Note that our data base includes both loosely (S n « 1 MeV) and well bound nuclei (S n ~ 2 — 8 MeV ) . For the loosely 
bound systems ( 14 B and 15 C ) R s = 0.8 ± 0.1 in agreement with the values deduced by Brown et al. [89] and Endcrs 
et al. [36] for 8 B and 9 C. 

B. Transverse momentum distributions: Carbon target 

The transverse momentum distributions (p x ) were calculated as described in Sections V and VI. In order to make 
a comparison with the experimentally measured distributions, the calculated distributions were filtered through a 
simulation which took into account the various broadening effects — straggling, the resolution of the tracking detectors 
and the spectrometer (Section II and Table I) — together with the finite angular acceptances of the spectrograph. The 
most significant effect was that arising from the acceptances, whereby the high momentum wings of the distributions 
were reduced in intensity 4 . Good agreement, as may be seen in Fig. 4, was found between the measured and calculated 
distributions. 

As mentioned in Section V-B, the transverse momenta are strongly influenced by the interaction with the target 
and the distributions arising from absorption and diffraction exhibit different lineshapes. This effect is explored in 
Fig. 14, where the longitudinal and transverse distributions have been calculated for s, p and d-wave single-particle 
states with binding energies of 1 MeV in an A=14 system. In the case of the longitudinal momenta both the s and 
(to a lesser extent) p-wave states result in relatively narrow distributions, whilst the d-wave state may be identified 
with a broad distribution exhibiting two symmetric peaks. 

In the case of the transverse momenta, the contribution arising from diffraction is systematically much broader than 
that from absorption and as such dominates the high momentum components of the total distribution. In addition, 
only the s-wave configuration gives rise to a relatively narrow distribution. The p and d-wave states are relatively 
broad with the former presenting a flat topped distribution with a small central dip. These features, combined 



3 The asymmetric nature of some of the distributions is discussed in the following section. 

4 As the emittance of the beam was relatively large (Section II), the angular acceptances of the spectrometer do not introduce 
a sharp cut-off in the transverse momentum distributions. 
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with the intrinsic structure of the projectile results in transverse distributions significantly different in form from the 
longitudinal distributions. In particular, the presence in a mixed configuration of a non negligible s-wave component 
will manifest itself in the transverse momentum distribution as a narrow feature superimposed on a much broader 
component. This is particulary well illustrated by the results for 23 and 24 ' 25 F (Fig. 4). Thus, whilst the lineshapes 
of the transverse momentum distributions may be more complex, they remain sensitive to the nature of the projectile 
ground state and may furnish spectroscopic information in a complementary manner. 

As noted in the preceding section a number of the longitudinal momentum distributions exhibit low momentum 
tails (Fig. 1). Such asymmetric distributions may arise from dissipative core-target collisions, such as observed in 
stable beam fragmentation [1,2,86] or, more likely in the case of weakly bound systems as a result of diffractive/elastic 
breakup [87]. Experimentally a correlation exists between the longitudinal and transverse momenta for events in the 
low momentum tail, as displayed in Fig. 15. Here the data for 15 C, where the yield is dominated by breakup to 
the 1 C ground state, was analysed so as to minimise any momentum shifts arising from core excited states. When 
events situated in the tail are selected the corresponding transverse momentum distribution is broad. In contrast, for 
events with p z greater than the mean momentum 5 the transverse momentum distribution is much narrower, with a 
width identical to that of the total distribution (Fig. 4). The low momentum events constituting the tail thus, on 
average, exhibit a much larger scattering angle. This result is consistent with the observation by Tostevin et al. [87] of 
asymmetric longitudinal momentum distributions at scattering angles away from zero degrees in the breakup of 15 C. 
Furthermore, the average momenta of the distributions were observed to be downshifted with increasing scattering 
angle. Such energy nonconservation effects cannot be described within the framework of the eikonal approximation 
employed here. As described in ref. [87], fully dynamical coupled discretised continuum channels calculations are 
capable of reproducing these effects, suggesting that the origin is diffractive/elastic breakup. A model independent 
confirmation could be furnished by fully exclusive measurements in which the beam velocity neutrons from breakup 
(a signature of diffractive dissociation) are measured in coincidence with the core fragments and deexcitation gamma 
rays. 

C. Longitudinal and transverse momentum distributions: Ta target 

As noted earlier, the longitudinal momentum distributions were measured for the breakup of 14 B, 15 ~ 17 C and 
17 ~ 19 N on Ta (Fig. 5), and distributions almost identical in width and form to those obtained for reactions on C were 
observed. Owing to the Coulomb deflection of the projectile and core fragment in the field of the target nucleus very 
broad transverse momentum distributions were encountered experimentally. As an example the distribution obtained 
for 15 C is displayed in Figure 16. It is clearly evident that the acceptances (p x ~ ±200 MeV/c) of the spectrograph 
were too limited to allow either transverse momentum distribution or the single-neutron removal cross section to be 
determined. 

An estimate of the effects of Coulomb orbital deflection on the dissociation of 15 C is shown in Fig. 16. For simplicity 
the transverse distribution has been calculated assuming pure Coulomb breakup (which is expected to dominate for 
the breakup of 15 C) and considering only the dominant s-wave component in the ground state wave function. The 
solid line in Fig. 16 corresponds to the assumption that the core fragment is deflected, following dissociation, along 
a classical Rutherford trajectory. In this case the angle of deflection depends only on the impact parameter b. The 
final transverse distribution was simulated assuming a distribution of impact parameters given by the calculations 
described in Section VI for b > b m i n = R core + Rtarget joined smoothly with a diffuse shape for b < b m i n mocked up by 
a Woods-Saxon form factor. The calculation in Fig. 16 shows that the broadening effect in the transverse momentum 
distribution is largely explained by core deflection in the target Coulomb field and suggests the breakdown of the 
straight line trajectory assumption at intermediate energies for heavy targets. 

We can now turn to the longitudinal core fragment momentum distributions on Ta target. Calculations including 
nuclear and Coulomb components compare well with the data as displayed in Fig. 5. As for the results obtained 
with the carbon target the theoretical predictions have been filtered through a Monte Carlo simulation (Section 
II) to take account the experimental effects. Details of the calculation are shown in Fig. 17 for 15 ' 17 C. In the 
case of 15 C, with a ground state dominated by an s-wave valence neutron and low S n , the nuclear and Coulomb 
distributions are almost identical, an example of the "numerical coincidence" pointed out by Hansen [17]. In the case 
of 17 C, which is dominated by a <i-wave valence neutron the Coulomb and nuclear distributions are not identical. The 
Coulomb interaction samples large impact parameters and selects small momentum components and the corresponding 



5 We note that such events should reflect most directly the intrinsic momentum of the removed neutron. 
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momentum distribution is narrow. For the cases detailed here, the laboratory grazing angle is about 1.4°, whilst the 
measured angular range is ±2°. Clearly a nuclear component must be present in these data and the calculations 
presented in Fig. 17 suggest that the absolute nuclear and Coulomb contributions predicted by the Glauber model 
coupled to first order perturbation theory for Coulomb dissociation seem to be realistic. 

D. Momentum distributions as a spectroscopic tool 

On the basis of the preceding comparisons, the reaction mechanism on the carbon target appears to be understood 
and, except for the low momentum tails, well described by a Glauber type approach within the eikonal approximation. 
When the nuclear structure is well known, as is the case for the nuclei closest to stability, the data are well reproduced 
by the model. In this section the manner by which the spin-parity assignments given in our earlier paper [39] to 
nuclei with poorly established ground state structure will be outlined. This will also be instructive in illustrating the 
sensitivity of the inclusive core fragment momentum distributions to the structure of the projectile. 

A comparison of the measured core fragment longitudinal momentum distributions from reactions on the carbon 
target and those calculated for the possible ground-state spin-parities is provided in Fig. 18. Here both the calculated 
and experimental distributions are displayed on an absolute scale (mb/[MeV/c]) without any normalisation, except 
for 23 O where no cross section could be extracted experimentally (Section III). In this case the calculated momentum 
distributions have been normalised so as to best reproduce the experimental distribution. In most of the cases the 
choice between the various JJ S is clear and the favoured spin-parity assignments, represented by the solid lines in 
Fig. 18, are listed in Table II 6 . Interestingly, in all the cases presented here, the favoured spin-parity assignments 
correspond to those suggested by the shell model calculations. 

In two cases ( 17 C and 23 0) the spin-parity assignments are not directly evident from inspection of Fig. 18. As 
it has been the object of recent attention [51,92,93] we will first turn our attention to 23 0. The form of the single- 
neutron removal longitudinal momentum distribution obtained here is well reproduced by both JJ s =l/2+ and 5/2+ 
assignments (similar results hold for the transverse momentum distribution). The former, however, leads to a predicted 
cross section of some 224 mb (Appendix B, Table IV) 7 , a factor of around 4 higher than for JJ s =5/2+. Based on 
the systematics of the cross sections obtained here for the Oxygen isotopes (Fig. 3), a 1/2+ assignment was favoured 
[39]. This conclusion is, contrary to the arguments made recently by the the authors of Ref. [51], confirmed by their 
measured cross section for single-neutron removal — 233±37 mb as compared to a value of 185 mb which we predict 
(Table IV), using the ground state structure given in Appendix B, at their beam energy of 72 MeV/nuclcon (a very 
similar result was also found by Brown et al. [92]). Further support for this ground state structure for 23 O may be 
found in the preliminary results for single-neutron removal at very high energy (936 MeV/nucleon) obtained using the 
FRS [53]. Not only can we reproduce well the measured longitudinal momentum distribution assuming a l/2 + ground 
state, but a cross section of 82 mb is predicted (Table IV ) as compared to the experimental value of 84±11 mb. 

In the case of 17 C the three possible ground state spin-parity assignments are shown in Fig. 18. It is clear that 
the l/2 + assignment grossly overestimates the cross section. The 3/2 + and 5/2 + assignments reproduce the data 
quite well, with the former providing a marginally better description of the central part of the distribution. As 
noted in Appendix B, a JJ s =3/2 + results in a large yield to the 16 C 2+ state in single-neutron removal from 17 C. 
The observation by Maddalena et al. [32] and Datta Pramanik et al. [94] of the corresponding 1.76 MeV gamma-ray 
transition thus confirms directly the 3/2 + assignment. 

Interestingly, the transverse momentum distribution measured here for 17 C also supports this assignment. This 
may be seen in Fig. 19 where the predictions (obtained with the same spectroscopic factors as above and without 
adjusting the parameters of the reaction calculation) for the three possible spin-parity assignments are compared to 
the measured distribution. This illustrates, as noted earlier, that the transverse momentum distributions can carry 
spectroscopic information complementary to that provided by the longitudinal momenta. 

In the light of the results described here, the evolution of the core fragment momentum distributions with T z may be 
understood, in particular through the competing contributions of the valence neutron ^2s!/ 2 and vld^/ 2 admixtures. 
In summary, following the crossing of N=8, the ground states of the N=9 - 14 B, 15 C - and N=10 isotones - 15 B, 16 C 
- are significantly influenced by the intruder z^2si/ 2 and ^2s 2 y 2 configurations, respectively. As the neutron number 



6 Owing to an error in the compilation of Ref. [45], the likely spin-parity assignments for 24 F were given as 1 + , 3 + rather than 
2 + , 3 + in our original paper [39]. 

7 The calculated cross section listed in Table 1 of our original paper [39] omitted the strong yield to the 3 + state predicted at 
around 4.8 MeV in 22 (Table IV). 
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increases so too does the contribution from vld§/2 configurations. This reaches a maximum at N=14 and then, as 
expected from the naive shell model, the the u2s\/2 orbital is occupied for N=15 and 16. Similar conclusions may 
also be drawn from the interaction cross section measurements of Ozawa et al., which exhibit enhancements for 23 O 
and 24 < 25 F [66,88]. 

IX. CONCLUSIONS 

An investigation of high-energy one-neutron removal reactions on 23 neutron-rich p-sd shell nuclei has been pre- 
sented. By studying isotopic chains extending from strongly bound near stable systems to weakly bound near dripline 
nuclei, the evolution of structure with isospin, as expressed by the core fragment observables (longitudinal and trans- 
verse momentum distributions and inclusive cross sections) has been explored. Experimentally, the measurements 
were carried out using a broad range high resolution, high acceptance spectrometer which permitted the data to be 
collected at a single magnetic field setting. Data were recorded using both Carbon and Tantalum targets in order to 
explore nuclear and Coulomb induced breakup. In the case of the Carbon target data, the large angular acceptances 
of the spectrometer were sufficient to encompass the full range of transverse momenta, thus permitting unambigu- 
ous measurements of the longitudinal and transverse momentum distributions and associated cross sections to be 
made. Owing to the large Coulomb deflection present in the reactions on the Tantalum target, the effective transverse 
momentum acceptances were very limited. As such only the longitudinal momentum distributions could be deduced. 

From the theoretical standpoint, an extended version of the Glauber model which incorporates effective NN in- 
teractions and second order noneikonal corrections to the JLM parameterisation of the optical potential has been 
developed. The treatment of Coulomb dissociation using first order perturbation theory has also been described. 
Particular emphasis has been devoted to retain the role played by the valence-neutron wave function via its Wigner 
transform in mapping the intrinsic momentum components onto the measured distributions. Despite a number of 
simplifying assumptions the model predictions agree very well with the experimental data, in particular those obtained 
for nuclei near stability with relatively well known structure. 

In the case of nuclear induced breakup the model suggests that for the longitudinal momentum distribution the 
reaction mechanism factorises in a manner such that only the low or surface momentum components in the wave 
function are selected. As a consequence only the asymptotic part of the wave function is probed. In the transverse 
momentum distributions this factorization does not occur and additional momenta arising from interactions with 
the target come into play. The principle drawback of the present approach, which is inherent to all Glauber-type 
models, is the neglect of energy and momentum conservation in describing diffractive/elastic breakup. The predicted 
momentum distributions are thus always symmetric and the low momentum tails observed here for some of the weakly 
bound nuclei cannot be reproduced. As noted in Section VIII-B, the description of such asymmetries requires the 
implementation of fully dynamical calculations (see, for example, [87]). 

Shell model spectroscopic factors calculated using the Warburton-Brown effective interaction formed the structural 
input for the calculations. The resulting momentum distributions and cross sections were found to be in very good 
agreement (except for the cross section for 22 F which was undcrprcdicted by some 30%) with the measurements. 
This agreement, especially for those nuclei with well established structure, suggests that the longitudinal momentum 
distributions and associated inclusive cross sections constitute a spectroscopic tool and ground state spin-parity 
assignments were proposed for 15 B, 17 C, 19 ~ 21 N, 21,23 , 23 _25 F. In addition to the dominance of the i/2si/2 intruder 
configuration in the N = 9 isotones, 14 B and 15 C, significant ^2s 2 ^ 2 admixtures were found to occur in the ground 

states of the neighbouring N = 10 nuclei 15 B and 16 C. Similarly, following the crossing the N=14 subshell, the 
occupation of the ^2s!/ 2 orbital is clearly observed for 23 , 24 < 25 F. 

The calculations of the transverse momentum distributions were also seen to agree well with the measurements. 
Thus, whilst being systematically somewhat broader than the longitudinal distribution the transverse distribution also 
carries structural information. Interestingly, due to the interplay of the projectile structure and reaction mechanism 
the transverse momenta were seen to often carry information in a complementary manner. In particular, the compe- 
tition between s and d— wave valence neutron configurations can exhibit itself directly in the transverse momentum 
distribution. Such complementary information may be of utility when conducting experiments with weak beams. 

Ultimately, the experimental determination of the core excited states populated in the reaction is required if detailed, 
unambiguous spectroscopic information is to be deduced. As seen elsewhere [28-34] this may be obtained for bound 
core states using large scale Nal or Germanium-detector arrays. As the neutron dripline is approached and the core 
itself becomes weakly bound, coincident neutron detection will also become necessary. 

Finally, in terms of perspectives, we conclude with some more general observations concerning the use of high- 
energy single-nucleon removal reactions as a probe of structure. As noted above, the reaction probes only the surface 
content of the projectile wavefunction. As such comparatively simple wavefunctions have been employed to describe 
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the valence nucleon. These wavefunctions, weighted by spectroscopic factors derived from large scale shell model 
calculations are, as described here, coupled to relatively sophisticated reaction model calculations. As is evident from 
the present work and has already been pointed out by others [27,35], remarkably good agreement has been achieved 
to date in describing the measurements. 

A few caveats should, however, be added. First, given the uncertainties inherent in the calculations, such as 
those described in Section VIII-A, uncertainties of order ±10% should be ascribed to the predicted cross sections 8 . 
Coupled with the experimental uncertainties — typically of a similar order — it would appear that deficiencies in our 
modelling of ^10-20% in cross section could easily be overlooked. High precision data obtained employing beams with 
very well established structure would, therefore, provide a means to help validate the accuracy to which the present 
approaches can be employed. A recent reanalysis of very high-energy inclusive measurements of single-nucleon removal 
from beams of 12 C and 16 [89] is an encouraging step in the this direction and dedicated experiments employing 
coincident gamma-ray detection are to be expected. 

Second, modelling employing "realistic" wavefunctions should be explored. In the case of single-nucleon transfer 
reactions it has long been known that despite their surface nature, the extraction of absolute spectroscopic factors 
can depend strongly on the description of the valence nucleon wavefunction [95]. In this context, it is instructive 
to recall a recent reanalysis of (d, 3 He) measurements by Kramer et al. [96]. In this study it was demonstrated that 
whilst only the tail of the bound-state wavefunction is sampled, it is very sensitive to the exact shape of the potential, 
thus introducing a significant model dependence in the calculated cross sections. In terms of weakly bound nuclei, the 
need to employ realistic wavefunctions was also found in the analysis of the p( n Be, 10 Be(2 + ))d reaction [97]. Similar 
effects must almost certainly be addressed in the analysis of high-energy, single-nucleon removal, and the levels to 
which they occur may provide a limit to the relatively simple analyses employed to date. 
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APPENDIX A 



In this appendix explicit analytical expressions for the Coulomb dissociation cross section Eq. (49) are given for 
s, p and (i-wave functions for E\, E 2 transitions and E\E 2 interference. The short hand notations FX X are defined 
below. Pi denote Legendre polynomials where the argument cos 6 has been omitted. The functions K ^ are defined in 
Section VI. 

FIei = k ^) + Ki{Oi 2 (Al) 

F2 E1 = -2K {0+Ki{Ol 2 (A2) 
Fl E2 = 3K Q (0 + K 2 (0 + Ki(0l 2 (2- P 2 ) 2 (A3) 

F2 E2 = 3^o(0-^2(0 + \Ki{i)l 2 (2 - d 2 ) 2 (A4) 

F3 E2 = 9tf (0 + 7^2(0 - 2^i(0l 2 (2 - I3 2 ) 2 (A5) 
Fl E1E2 = 2K (Z)+K 1 (Z) 1 2 (2- (3 2 ) (A6) 
F2 E1E2 =3K (0- K 1 (Ch 2 (2-/3 2 ) (A7) 
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L=2 
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APPENDIX B 



In the following, listings are provided of the calculated spectroscopic factors (C 2 S) and cross sections (a(I^)) to 
the core excited states (E^ x , If?) populated in single-neutron removal from the projectile nucleus ( A Z, J 7r ) by the 
Carbon target. The contributions arising from absorption (a a b s ) and diffractive dissociation (cdiff) are detailed 
(note: the latter includes the very small contributions arising from Coulomb breakup) and the total inclusive cross 
section (<T^^ l t ' 6er ) to all core states is given. Projectile JJ S and core excited state energies taken from the shell model 
predictions are marked in parentheses. Only the results for the preferred ground structure are given here: the listings 
for other J 71 ' may be found in rcf. [45]. 
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22.9 






0.953 2+ 


lPl/2 


1.17 


22.6 


16.7 


39.3 












iJLauber 
"-In — 


62 mb 


14 B 


2" 


g.s. 3/2" 


ld 5 / 2 


0.31 


9.6 


8.7 


18.3 








2si/ 2 


0.64 


57.0 


64.7 


121.7 






3.483 3/2+ lp 1/2 


0.41 


8.5 


7.1 


15.6 






3.68 5/2+ 


lPl/2 


0.8 


16.2 


13.3 


29.5 














:ber =1 


85 mb 


i0 B (3/2") 


g.s. 2- 


ld 5 /2 


0.28 


5.8 


5.3 


11.1 








2si/ 2 


0.48 


20.3 


22.3 


42.6 






(0.89) 1" 


2si/ 2 


0.27 


9.4 


9.9 


19.3 






(0.73) 1+ 


lPl/2 


0.58 


11.7 


11.3 


23 






(0.96) 3" 


ld 5 / 2 


0.47 


8.7 


7.6 


16.3 


(T (Jlauber = 1U mb 



25 





E^ x 


11 


nlj 


C' 2 S 


&abs &diff &(Ir) 




[McV] 








[mb] [mb] [mb] 


i4 C 0+ 


g.s. 


1/2" 


lPl/2 


1.67 


26.4 15.9 42.3 




3.089 


1/2+ 








3.684 


3/2" 


1P3/2 


2.05 


29.7 16.9 46.6 


a Ulauber =89 mb 


i0 C 1/2 + 


g.s. 


0+ 


2Si/2 


0.83 


62.1 62.1 124.2 




6.094 


1" 


1P3/2 


0.16 


2.8 1.9 4.7 








lPl/2 


1.03 


16.3 10.9 27.2 




6.903 


o- 


lPl/2 


0.46 


6.9 4.6 11.5 


(T G(at l 6er = 168 mb 


iy C 0+ 


g.s. 


1/2+ 


2si/ 2 


0.6 


19.3 17.5 36.8 




0.740 


5/2+ 


ld 5 / 2 


1.23 


21.9 16.2 38.1 




iY C (3/2+) 


g.s. 


0+ 


ld 3 / 2 


0.035 


0.9 0.8 1.7 




1.762 


2+ 


ld 5 /2 


1.41 


29.3 25.5 54.8 








2si/ 2 


0.16 


6.9 7.2 14.1 




4.1 


2,3,4+ ld 5/2 


0.76 


12.5 10.0 22.5 








2si/ 2 


0.22 


6.1 5.7 11.8 


(T U(a«ber = 105 mb 


i8 C + 


g.s. 


3/2+ 


ld 3 / 2 


0.1 


1.4 1.2 2.6 




(0.04) 


5/2+ 


ld 5/2 


2.8 


43.3 38.0 81.3 




(0.3) 


1/2+ 


2si/ 2 


0.65 


17.2 17.7 34.9 


a (Jlauber^ U Q mb 



A z r 


E c ex 


r 


nlj 


C 2 S 


Cafes a diSf 






[MeV] 








[mb] [mb] 


[mb] 


iV N 1/2- 


g.s. 


2" 


ld 5 /2 


0.59 


11.1 7.6 


18.7 




0.120 


0" 


2si/ 2 


0.12 


3.4 2.7 


6.1 




0.298 


3" 


ld 5 /2 


0.784 


14.5 9.8 


24.3 




0.397 


1" 


ld 5 / 2 


0.36 


9.9 7.8 


17.7 












—(Jtauber 

a -ln - 


67 mb 


i«N 1- 


g.s. 


1/2" 












1.374 


3/2- 


ld 5/2 


0.65 


13.4 10.1 


23.5 








2si/ 2 


0.195 


6.7 5.9 


12.6 




1.850 


1/2+ 










1.907 


5/2- 


ld 5 /2 


0.89 


15.4 11.4 


26.8 




3.129 


7/2- 


ld 5 /2 


0.42 


7.6 5.4 


13 




3.2 


3/2- 


2si/ 2 


0.15 


4.1 3.4 


7.5 








ld 3 /2 


0.29 


4.6 3.2 


7.8 












—Glauber 
"-In — 


91 mb 


iy N (1/2-) 


g.s. 


1- 


ld 3 /2 


0.02 


0.3 0.24 


0.54 








2si/ 2 


0.005 


0.14 0.16 


0.3 




0.115 


(2-) lds/2 


1.26 


21.6 16.9 


38.5 




0.747 


(3") 




1.71 


24.9 18.7 


43.6 




(0.936) (1-) 


2S1/2 


0.35 


8.6 7.7 


16.3 












—Glauber 
"-In — 


99 mb 


a, N (2-) 


g.s. 


1/2" 


ld 5 /2 


0.36 


8.5 7.9 


16.4 




(1.68) 


3/2- 


ld 5 /2 


0.66 


12.4 10.7 


23.1 




(2.17) 


5/2- 


ld 5 / 2 


0.38 


6.8 5.8 


12.6 




(3.9) 


7/2" 


ld 5 / 2 


1.73 


26.8 21.6 


48.4 


(T (Jlauber = 1Q1 mb 


ai N (1/2-) 


g.s. 


2- 


ld 5 /2 


1.744 


28.3 25.1 


53.4 




(0.6) 


3- 


ld 5 /2 


2.61 


40.1 35.1 


75.2 




(0.74) 


1- 


2si/ 2 


0.45 


11.1 11.1 


22.2 



r Ef x T* nlj C Z S a abs a dtff a(I*) 
[MeV] [mb] [mb] [mb] 

iy 5/2+ g^ 0+ ld 5/2 0.685 14.1 9.8 23J~ 

1.982 2+ ld 5/2 0.48 8.4 5.6 14.0 

ld 3/2 0.019 0.3 0.2 0.5 

2s 1/2 0.009 0.25 0.2 0.45 

3.555 4+ ld 5/2 1.24 19.7 12.6 32.3 

3.92 2+ ld 5/2 0.22 3.4 2.2 5.6 

2si /2 0.06 1.34 0.96 2.3 

5.25 2+ ld 3/2 0.023 0.3 0.18 0.48 

ld 5/2 0.016 0.2 0.18 0.38 

5.38 3+ 2si /2 0.11 2.2 1.5 3.7 

cr G ^ tlbe7 =84 mb 

20 0+" g.s. 5/2+ ld 5/2 3.43 51.3 34.5 85.8 
0.096 3/2+ ld 3/2 0.05 0.66 0.44 1.1 
1.471 1/2+ 2si /2 0.28 5.3 4.0 9.3 

cr G ^ Mber =96 mb 

(5/2+) g.s. 0+ ld 5/2 0.345 6.7 5^4 127T 

1.67 2+ ld 5/2 1.3 21.8 16.4 38.2 

2si /2 0.004 0.1 0.1 0.2 

3.57 4+ ld 5/2 2.59 38.0 27.3 65.3 

4.072 2+ ld 5/2 0.09 1.3 0.9 2.2 

2s 1/2 0.05 1.0 0.9 1.9 

4.456 0+ 
4.85 4+ 

5.23 2+ ld 5/2 0.12 1.6 1.1 2.7 

<7 Gl 1 a r ^ ber ^123 mb 

0+" g.s. 5/2+ ld 5/2 5.22 74.9 56.9 131.8 

1.33 1/2+ 2si /2 0.23 4.3 3.7 8.0 

2.2 3/2+ ld 3/2 0.03 0.33 0.27 0.6 

3.08 5/2+ ld 3/2 0.14 1.7 1.2 2.9 

o-^ 6er = 143 mb 



^ 3 (1/2+) g.s. 


0+ 


2si/ 2 


0.8 


31.7 


33.8 


65.5 


3.38 


2+ 


ld 3 / 2 


0.053 


0.7 


0.5 


1.2 






ld 5 /2 


2.1 


30.5 


24.5 


55.0 


(4.62) 


0+ 


2si/ 2 


0.11 


2.2 


1.9 


4.1 


(4.83) 


3+ 


ld 5 /2 


3.08 


40.5 


31.7 


72.2 


(6.5) 


2+ 


ld 5 /2 


0.24 


2.9 


2.1 


5.0 


(6.64) 


o- 


lPl/2 


0.36 


3.4 


2.6 


6.0 


(6.9) 


1" 


lPl/2 


0.94 


8.6 


6.5 


15.1 



Gfg M ber =2 24^ b - 
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A z 


J 11 


E c 

ex 


r 


?li 7 


C' 2 S 


'-' a,os 










\MeV] 








[mbl 

L J 


[mb 

L J 


[mbl 

L J 




4+ 


g.s. 


5/9+ 

O 1 Z 


IQ5/2 


n ^r 




4 

U . L ± 


IO.O 






9799 


1 /2+ 

1/ z 














1. / o 


■} /9 + 
o/ Z 


1 rl , 


n 9/t 
u.z^ 


0.0 


9 A 


o.u 






C1 9"l 


Q/9+ 
y/ z 


ld 5 / 2 


n or 

u.yu 


10. U 


Q 1 


99 7 






f3 56) 

I O . UU J 


3/2+ 

O / Z 


1^5/2 


n nQ 

U . 


1 1 Q 


U. O -L 


9 






f3 64) 


7/2+ 
//z 


1^5/2 


n ns7 


1 1 5 


75 

U. 1 O 


1 9 












U.U10 


U.o 


U.z 


u.o 








Q/9+ 


lds/2 


u .00 


7 4 


4 & 


1 9 9 






(4.U2J 


1/2 














I'd 4^ 
l+- 40 7 


7/9+ 
'/ z 


J-Q-5/2 


0^ 

U . UO 


4 

U.I 


9 
u . z 


u .u 










9s-, m 
zb l/2 


09 

U . UZi 


u.o 


9 

U.Zi 


^ 
u . 








1 1 /9+ 
± 1/ z 


1 rl , 

id 5/2 


U.04 


0.0 


/I 1 


1U.0 










ld3/2 


04 

VJ . VJ4: 


4 


^ 
u.o 


7 

U . f 








13/2+ 


1 rl r- In 

±u 5/2 


0.52 


6.3 


3.9 


10.2 






(5 31) 


7/2+ 


lflr /o 
ila 5/2 


0.07 


0.8 


0.5 


1.3 










^ s l/2 


0.03 


0.5 


0.3 


0.8 






(5.5) 


7/2+ 


lcU /9 
J -^ l 0/2 


0.012 


0.2 


0.1 


0.3 






(6.95) 


13/2+ 


1 rl f 
lu 5/2 


0.075 


0.8 


0.5 


1.3 














—(Jiauber 
"-In — 


87 mb 


' za F (5/2+) 


g.s. 


4+ 


1 rl , 


1 9 

l.Z 


15.9 


11.1 


97 n 








o 


1 rl , 


n 7R 


10.0 


6.8 


IRS 
10.0 










9s-, ;o 


06 

u. uu 


1.1 


0.9 


9 

Zi .U 






en 7) 


2+ 


ld.5/2 


64 

U. U4: 


8.1 


5.6 


1 3 7 










9s, ,„ 


O 1 } 


0.63 


0.49 


1 1 9 

X . J.Z- 








5+ 

U 


lds/2 


1 4 

_L .4: 


17.0 


11.4 


98 4 

ZiO .4: 






(\ 6) 


1 + 


- LU 5/2 


0.06 


0.7 


0.5 


1.2 






(1.65) 


3+ 


2si/ 2 


0.07 


1.2 


0.8 


2.0 






(1.67) 


2+ 


ld 5 / 2 


0.31 


3.7 


2.5 


6.2 






(2.3) 


1+ 


ld 5 / 2 


0.27 


3.1 


2.1 


5.2 






(3.5) 


5+ 


ld 5 / 2 


0.11 


1.2 


0.8 


2.0 


a Ulauber = 10 Q mb 



^ 4 F (3+) g.s. 


5/2+ 


ld 5 / 2 


0.09 


1.6 


1.3 


2.9 






2si/ 2 


0.74 


22.3 


21.1 


43.4 


(1.8) 


1/2+ 


ld 5 /2 


0.073 


1.1 


0.8 


1.9 


(2.9) 


7/2+ 


ld 5 /2 


0.44 


6.0 


4.4 


10.4 


(3.2) 


5/2+ 


ld 5 /2 


0.37 


4.9 


3.6 


8.5 


(3.7) 


9/2+ 


ld 5 /2 


0.96 


12.3 


9.0 


21.3 


(4.2) 


7/2+ 


ld 5 / 2 


0.38 


4.7 


3.4 


8.1 


(4.4) 


3/2+ 


ld 5 /2 


0.2 


2.5 


1.7 


4.2 


(4.61) 


5/2+ 


!d 5 /2 


0.09 


1.1 


0.8 


1.9 






2si/ 2 


0.04 


0.7 


0.5 


1.2 


(4.65) 


9/2+ 


ld 5 /2 


0.2 


2.4 


1.7 


4.1 


(4.78) 


5/2+ 


2Sl/2 


0.052 


0.9 


0.7 


1.6 


(5.6) 


5/2+ 


ld 5 /2 


0.12 


1.4 


0.9 


2.3 


(6.77) 


3/2+ 


ld 5 /2 


0.16 


1.7 


1.2 


2.9 



^Glauber mb 
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A z r 


ex 


H nlj 


C' 2 S 






trim) 

\ C I 




[MeVl 

L J 






[mb] 

L J 


[mb 

L J 


[mb] 

L J 


' zi> V (^l2+\ 


g.s. 


3+ 9«, 


n 89 


91 


91 


49 






ia 5/2 


U.Uo 


1 ^ 

l.o 


i n 

l.U 


9 9 
Z.o 




rn 1 

(U.1J 


9+ 9«, /„ 
z zb l/2 




1U.U 


1 f{ ft 


^9 




(<J.O) 


1 iu 3/2 


flfi 
u.uu 


7 




1 s 




(2 2) 


4+ 1Hc/n 

^ ±u 5/2 


1.0 


13 


10.0 


23 




(2.5) 


3+ ld 5/2 


0.7 


8.4 


6.3 


14.7 




(2.8) 


1+ ld 5/2 


0.18 


2.1 


1.6 


3.7 




(2.9) 


1+ ld 5/2 


0.1 


1.2 


0.9 


2.1 




(3.5) 


5+ ld 5/2 


1.6 


18.3 


13.5 


31.8 


a (Jlauber = lM mb 
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TABLE I. Examples of the contributions of different experimental broadening effects on the widths (FWHM) of the core 
fragment momentum distributions (see text). 





A Z 




Lorentz 


target 




spectrometer 




broadening 






intrinsic 








=> 


measured 


[%] 


FWHM pz 




56.3 


59 




63 




63.6 


11 


(C target) 


19 N 


166 


175.5 




187.6 




188.4 


11.6 


FWHM pz 




57 


60 




61.8 




62.2 


8.4 


(Ta target) 




176 


186 




193 




194 


9.3 






A Z 




target 


tracking 




spectrometer 




broadening 






intrinsic 










measured 


[%] 


FWHMpj 




75 


77 




76.8 




79 


5 


(C target) 


19 N 


226 


227 




228 




229 


1.3 
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TABLE II. Summary of the results obtained with the Carbon target. Where available the results of other experiments are 
also listed. 



Ay 


energy 


r vv jrLivip 2 


r vv nivip z 


r vvnivipj. 




f-ln 


_.Ulauber 
G -ln 


JIT 




[Me V / nucleon] 


[Me V / c| 


[Me V / cj 


[Mev / Cj 


[Me V / c| 


[mbj 


[mbj 




1!2 B 


64 


158 ± 3.5 


142 ± 3.5 


175 ± 3 


173 ± 3 


81 ± 5 


91 


1 + 


13g 


54 


150 ± 7 


135 ± 7 


178 ± 2 


176 ± 2 


59 ± 4 


62 


3/2" 


14-r> 

rJ 


00 

86 
59 


CO £? _1_ A g. 

OO.O ± U.O 


r/; rr _r_ n rr 
00.0 ± 0.0 

57 =b 2 a 
55 ± 2 6 


<y ± 2 


(0 ± 2 


iro I i rr 

loo ± lo 
48 ± 5 a 
176 =b 16 b 


1 Q g. 

loo 


o — 

Z 


15 R 


A Q 

4o 


CO 1 O E 
OZ zt Z.O 


16 ± 2.0 


84 ± y 


on I n 

oU ± y 


i no i i ■ i 
lUo it lo 


i in 

llz 


o fry— C 
6/2 


14 C 


67 


200 ± 5 


180 ± 5 


222 ± 3 


220 ± 3 


65 =b 4 


89 


0+ 


15 C 


62 

54 
85 


71 ± 0.7 


63.5 ± 0.7 
67 ± 3 a 


81 ± 1 


86 ± 1 


159 =b 15 
137 ± 16 d 

33 ± 3 a ' q 


168 


1/2+ 


i6ri 


r r 

00 

62 
oo 


121 ± 2 


i no _i_ o 
108 ± 2 

yu m y 


14o ± O 


1 A C\ _l_ Q 

14U ± 6 


00 ± 

77 =b 9 e 

+ r 
00 -10 


(0 





17p 
L/ 


/IQ 
4y 

o4 
02 
yo.o 

Q0.4 
yu4 


120 in 


111 4- *3 
ill ± o 

1 /I E _1_ E a 

140 ± 

y4 it ly 

1/M 4- fi 9 

141 it 


ioy zt y 


i aa _i_ a 
100 zt y 


OA _[_ Q 

o4 zt y 

Oft. -\- Q a ' ( ? 
ZO Zt O 

1 1 g, _l_ 1/1 e 
110 zt 14 

/II _l_ /I /f9 
41 ±4 

1 9Q-I-99 h 
izyztzz 


1UO 


o /o+ c,a,e,o,t 
of Z 


i8ri 


40 
86.2 


140 It 


1 Oft 4- ^ 
120 ± 

1 1 A _J_ 1 O / 

11U ± 12 


1 CQ J_ 1 ft 

ioy zt id 


1 _|_ i 
100 zrz 10 


11C _|_ 1 Q 

110 zt lo 

or _l_ o J?*? 
,30 ± Z 


1 1 Q 

iiy 


n+ 
u 


17 AT 

JN 


65 


158 ± 4 


141 ± 4 


217 zb 4 


214 zt 4 


55 zb 5 


67 


1/2 


JN 


59 


188 ± 3 


168 ± 3 


219 =b 3 


216 zb 3 


109 zb 11 


91 


1~ 


IN 


53 


199 ± 3 


177 ± 3 


229 ziz 5 


226 zb 5 


86 ± 9 


99 


1/2 ' M 


20 m 
IN 


48 


184 ± 4 


162 ± 4 


220 ± lb 


217 dz 16 


98 zb 13 


101 


2 


21 N 
IN 


40 


1 i Olt 1 


1 /IQ 4- 7 

i4y it ( 






14U zt 44 


101 


1/2 


(J 


68 


214 ± 8 


190 ± 8 


253 ± 13 


250 ± 13 


104 d= 12 


O A 

84 


5/2 


20 r> 
(J 


62 


247 ± 5 


219 ± 5 


254 ± 7 


251 it 7 


112 zb 11 


96 





(J 


56 


237 ± 6 


210 ± 6 


246 ± 7 


243 it 7 


134 d= 14 


123 


5/2 ,H 


22q 


51 


235 ± 4 


206 ± 4 


240 ± 16 


237 ± 16 


120 zb 14 


143 


o+ 


23q 


47 


135 ± 9 


114 ± 9 


162 ± 28.5 


157 ± 28.5 


TL 


224 3 


1/2+ c 




72 


94 ± 12 j 


175 ± 14 3 






233 ± 37 3 


185 s 






936 




123 ± 10 p 






85 =b 15 p 


82 s 






64 


212 ± 14 


185 ± 14 


278 ± 28 


274 ± 28 


121 ± 16 


87 


4+ 


23 p 


59 


267 ± 4 


235 ± 4 


236 ± 10 


232 ± 10 


114 ± 12 


106 


5/2+ c ' fc ' ! 


24 p 


54 


151 ± 4 


129 ± 4 


203 ± 18 


198 =b 18 


124 ± 16 


115 




25p 


50 


128 ± 8 


106 ± 8 


173 ± 45 


168 =b 45 


173 =b 46 


154 


5/2+ c 



a ref. [49] (Be target), b ref. [30] (Be target), 

c assignment from present experiment, d ref. [87] e ref. [32] (Be target), / ref. [48] (Be target), 
g ref. [38], (C target), h ref. [52] (C target), 

i ref. [81], j ref. [51] (C target), k ref. [82], / ref. [83], m ref. [84], 

n no beam intensity normalization available, o ref. [85], p ref. [53] (C target), 

q affected by limited transverse momentum acceptances, r ref. [50], s see Table IV, t ref. [94] 
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TABLE III. Summary of the results obtained with the Tantalum target. Where available the results of other experiments 
are also listed. 



A z 


Encrffv 


FWrlMp" 6 


FWHMp^ 


(J 


Glauber 
® — In 




[MeV /nucleon] 


[MeV/c] 


[MeV/c] 


[mbl 


[mbl 




50 
59 
86 


62 ± 2 


57 ± 2 
59 ± 3 b 
48 ± 3 a 


ti 

638 ± 45 6 
157 a ' e 


864 


15 C 


62 
85 


69 ± 0.5 


63 ± 0.5 
67 ± 1 a 


_ 

75 a ' e 


978 


16 C 


55 


106 ± 3 


97 ± 3 


_ d 


193 


17 c 


49 
61 
62 


131 ± 7 


121 ± 7 
-110 y 


_ d 

350 ± 67 c 


280 


iv N 


65 


147 ± 5 


134 ± 5 


_ d 


173 


18 N 


59 


176 ± 5 


159 ± 5 


_ d 


238 


19 N 


53 


194 ± 11 


176 ± 11 


d 


216 


a ref. [49] (Ta target), b ref. [30] (Au target), 


c ref. [79] (Au 


target), 





d no reliable cross section could be estimated owing to very broad transverse momentum distributions (see text), 
e affected by restricted transverse momentum acceptances, f ref. [80] ( 238 U target) 



TABLE IV. Predicted cross sections for one-neutron removal from 23 O (Jg S = l/2 + ) by a Carbon target at 47 (present 
experiment), 72" and 900 6 MeV/nucleon. The core excited states and spectroscopic factors are those listed in Appendix B. 



E% [MeV] 


II 


47 


Eiab [MeV/nucleon] 
72 


900 


g.s. 
3.38 

4.62 
4.83 
6.5 
6.64 
6.9 


0+ 

2+ 

0+ 
3+ 
2+ 
0" 
1" 


65.5 
1.2 
55 
4.1 
72.2 
5.0 
6.0 
15.1 


53.0 
0.99 
45.4 
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FIG. 1. Comparison of the core fragment longitudinal momentum (p z ) distributions obtained using a Carbon target and the 
Glauber model calculations (solid line). 
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FIG. 2. Comparison of the widths (FWHM) of the longitudinal (filled circles) and transverse (open circles) core fragment 
momentum distributions for reactions on the Carbon target. 
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FIG. 3. Experimental one-neutron removal cross sections (filled circles) compared to the results of the Glauber calculations 
- open circle: total cros section, open triangle: absorption and open diamond: diffraction (see text for details). The points are 
connected by lines to guide the eye. 
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FIG. 4. Comparison of the core fragment transverse momentum (p x ) distributions obtained using a Carbon target and the 
Glauber model calculations (solid line). 
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FIG. 6. Coordinate system used in the Glauber model calculations. 
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FIG. 7. RMS matter radii extracted from Hartree-Fock + BCS calculations of single particle densities (black line), compared 
with experimental values from Liatard et al. [62] (filled circles), Ozawa et al. [66] (open circles) and Tanihata et al. [67] (open 
diamonds). The relativistic mean field calculations of Ren et al. [63-65] are indicated by the dashed lines. 
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FIG. 8. S-matrix, absorption profiles and transmission coefficients for n- 12 C (left panels) and 14 C- 12 C (right panels) inter- 
actions at 30 MeV/nuclcon obtained from the JLM folding model calculations. Results using the lowest order of the eikonal 
approximation are shown as a dotted line, those using the first order corrections by the dot-dashed line, and those using second 
order corrections by the solid line. For 14 C+ 12 C the calculations converge with the use of first order corrections. 
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FIG. 9. The effect of noneikonal corrections on the distortion kernels (see the text for definitions) as a function of impact 
parameter s, for neutron absorption and diffraction, core ( 14 C) absorption and for absorption of both the core and neutron 
by 12 C at 30 MeV/nuclcon. Dotted lines: lowest order eikonal approximation. Dashed-dotted lines: first order noneikonal 
corrections included. Solid line: second order corrections included. 
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FIG. 10. Elastic (dotted), reaction (dash-dotted) and total (solid line) neutron cross sections for reaction on four targets as 
a function of incident energy. The calculations were performed using JLM microscopic potentials in the eikonal approximation. 
Non eikonal corrections up to second order were included and the effects are shown by thin continuous lines (the convention is 
as per the previous two figures) for Be, C and Si targets. For the Ta target, the eikonal series does not converge in the range 
of energies shown. The experimental data for the total, elastic and reaction cross sections (open and closed circles and open 
squares) were taken from refs. [69,70] 
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FIG. 11. Selected examples (see text) of the core fragment longitudinal momentum distributions obtained using the Carbon 
target showing the contributions from the different reaction mechanisms. The calculated distributions (thick solid lines) include 
the absorption (thin solid lines), diffraction (dashed) and Coulomb (dash-dotted) components. 
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FIG. 12. Calculated absorption (solid lines) and diffraction (dashed lines) cross sections versus binding energy for s and 
d-wave states of an A— 17 system at 50 MeV/nucleon on a Carbon target using the JLM interaction and various orders of 
eikonal theory. 
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FIG. 13. Ratio of experimental to theoretical cross sections (R s ) as a function of projectile mass number for the data obtained 
with the Carbon target. The dotted line indicates the mean value and the shadded band the la variance. 
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FIG. 14. Test case calculations of the longitudinal (k z ) and transverse (k x ) core fragment momentum distributions for 
one- neutron removal by a Carbon target at 50 MeV/nuclcon assuming an s, p or d-wave state in an A=14 system with S„= 1 
MeV. The total (solid line), absorption (dashed), diffraction (dash-dotted) and Coulomb (dotted) components are indicated. 
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FIG. 15. Correlations between the longitudinal and transverse 1 C core fragment momenta for one-neutron removal from 
C by the Carbon target. 
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FIG. 16. Core fragment transverse momentum distributions for one-neutron removal from 15 C by the Tantalum target (filled 
circles). Dashed line: calculated transverse momentum distribution for Coulomb induced breakup. Solid line: same calculation 
after convolution with the orbital deflection in the Coulomb field of the target (see text for details) . 



47 




FIG. 17. Examples of the core fragment longitudinal momentum distributions for reactions on the Ta target. The filled 
circles represent the data. To aid in the comparison the calculated distributions for the total (solid lines) , nuclear (dotted) and 
Coulomb dissociation (dashed line) are normalised to the data. 
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FIG. 18. Comparison of the measured core fragment longitudinal momentum distributions (from reactions on the Carbon 
target) with those predicted for different projectile ground state spin-parity assignments. The favoured assignements are 
displayed as solid lines (see text). 
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FIG. 19. Comparison of the measured C core fragment transverse momentum distribution for reactions of C on the 
Carbon target with the distributions calculated for the three possible ground state spin-parity assignments. 
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